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Abstract 

This paper presents a six-band k • p theory for wurtzite semiconductor nanostructures with cylin- 
drical symmetry. Our work extends the formulation of Vahala and Sercel [Physical Review Letters 
65, 239 (1990)] to the Rashba-Sheka-Pikus Hamiltonian for wurtzite semiconductors, without the 
need for the axial approximation. Results comparing our formulation for studying the electronic 
structure of wurzite quantum dots with the conventional formulation are given. 
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I. INTRODUCTION 



The III-V nitride semiconductors, with wurtzite (WZ) crystal structure, have received a 
great deal of attention in recent years. In the 1950 's there were numerous optical studies of 
bulk WZ semiconductors (see Ref. Qfor a recent review). However, the detailed band struc- 
ture of these materials was not studied until the discovery (in 1993) of blue light emission 
of WZ GaN on sapphire.— Currently the interest in WZ materials has shifted to nanostruc- 
tures such as CdSe quantum rods,^»^ ZnS nanowires,— ZnO nanorods^ and nanowires,- 8 A1N 
nanorods^ and GaN nanowires.— On the theory side, the k ■ p theory of WZ bulk materials 
was developed by Rashba and Piku o 11 ' 12 and later applied to heterostructures by a number 
of authors . 13 i 14 Mireles and Ulloa^ have applied the theory to heterostructures using the 
envelope function formalism of Burlj^ and Foreman-^ As far as we are aware, the model 
has only been applied to quantum wellsi 8 - and pyramidal quantum dots.— In 1990, Sercel 
and Vahala (SV) presented a new formulation of the multiband envelope function theory 
and applied it to spherical quantum dots and cylindrical quantum wires of zincblende (ZB) 
materials . 20 ! 21 For ZB, the SV formulation was only possible provided the axial approxima- 
tion (72 = 73) was made in the Hamiltonian. The theory, within the Luttinger-Kohn (LK) 
framework, has subsequently been applied to model quantum rods,— quantum rings,— and 
quantum dots.— 

The current interest in WZ nanowires^ motivates the need for a model similar to the one 
introduced by Sercel and Vahala for ZB, for increasing the physical understanding of the 
k • p theory and for improving the efficiency of numerical computations for device applica- 
tions. In this Letter, we propose a formulation of the Rashba-Sheka-Pikus Hamiltonian with 
cylindrical symmetry. We have reformulated the k • p Hamiltonian in terms of the Sercel- 
Vahala (SV) representation for problems with axial symmetry. The SV representation is 
useful because it reduces the 3D problem to a 2D one when cylindrical polar coordinates are 
used. Contrary to the work of Sercel and Vahala, where the axial approximation had to be 
introduced for ZB materials, we show that the formulation is exact for WZ materials, i.e., 
no axial approximation was needed. In addition to the fundamental interest in the SV repre- 
sentation of the WZ Hamiltonian, the latter also helps make for more efficient computation 
for problems with axial symmetry. This includes free-standing and embedded cylindrical 
nanowires, modulated nanowires, quantum rods, spheroidal and spherical quantum dots. 
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II. HAMILTON! AN 



The six-band Hamiltonian for a WZ semiconductor heterostructure, in the \X f), \ Y f 
), \Z t), \X |), \Y |), |Z 1) basis states of the Kane model i a 12 i 14 i 15 i 25 
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troduced by Mireles and UUoa*^; the Aj's are spin-orbit parameters.— Next, one can use 
the \u\), \u 2 ), \u 3 ), \ui), {us), \uq) basis states defined in Refill. The advantage of using this 
basis is that now the zone-center Bloch functions are described by an angular momentum J. 
The envelope part of the total wave function behaves like the spherical harmonics, with an 
angular momentum L. The periodic Bloch space and the slowly varying envelope space are 
coupled by the k ■ p interaction. For problems with cylindrical symmetry, the projection of 
the total angular momentum F z = L z + J z is a good quantum number. 

Next, we express the k ■ p Hamiltonian in terms of cylindrical polar coordinates p, 0, z. 
The total wave function is then written as 
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There is a double degeneracy with respect to the sign of F z due to inversion and time-reversal 
symmetry. The new matrix elements are given by 

</<(r)|%i(r)> 

= (g i (p,z)\e^ F '- J «»HeW- J '>»\g j (p,z)) 

= {g i (p,z)\H'\g j (p,z)). (5) 

These matrix elements are zero unless F z = F' z . 

It is found that all <fi dependence goes away following the SV transformation, without 
making an axial approximation as for ZB. The validity of this result is due to the axial 
symmetry already found to be true for the bulk dispersion relation.— This can also be seen 
from a group theoretic point of view by noting that the group of the wave vector is the same 
for all wave vectors in the plane perpendicular to the c axis. Finally, the WZ Hamiltonian 
in cylindrical coordinates is 
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where J n = J zn . Note that in general ^ S^. That the matrix representation is nonhermi- 
tian is due to the chosen representation and to the fact that the Hilbert space is now divided 
into F z subspaces. An analogous result was obtained previously for the SV representation 
of zincblende quantum rings.— 



III. NUMERICAL RESULTS AND DISCUSSION 

A free standing WZ GaN cylindrical quantum dot, with a radius of 50 A and a height of 
100 A, with the axis along the crystallographic c-axis, has been studied using both the 3D 
and 2D methods. The parameters for GaN have been taken from the literature— The 3D 
Hamiltonian (Q) and the 2D Hamiltonian (jHJ) have been implemented using FEMLAB; this 
is a software using the finite element method. Different meshes have been generated, and 
the convergence of the eigenvalues has been studied. Due to the reduced dimensionality, 
the eigenvalue problem in 2D can be solved much faster, in minutes instead of hours, and 
with less memory requirements than the eigenvalue problem in 3D. After only one mesh 
refinement (2040 elements) the 2D eigenvalues have converged. However, even after three 
mesh refinements (3967 elements) the 3D eigenvalues are not fully converged, as seen from 
Fig.ffl 

A free standing WZ CdSe cylindrical quantum dot, with a radius of 50 A and a height of 
40-150 A, has also been investigated. The parameters for CdSe have been taken from the 
literature . 13 i 28 By varying the height we can study shape effects on semiconductor nanocrys- 
tals. From Figures we notice that there are crossings between states with F z = 1/2 
(solid lines) and F z = 3/2 (dotted lines). The inclusion of the linear term in the Hamilto- 
nian dramatically changes the energy band structure, and results in anticrossings between 
energy bands^. 
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In CdSe nanocrystals with a radius of 15 A and a gelcap shape a crossing between the 
lowest two valence sublevels appears at an aspect ratio of 1.25.— There are also some ex- 
perimental studies of the transition from 3D to 2D confinement in CdSe quantum rods.^i 
Theoretical work on the electronic properties of the transition reported so far include em- 
pirical pseudopotential calculations on a rather artificial shaped and k ■ p calculations for 
cylindrical ZB structures*^ and for spheroidal WZ structures.-^ We hope that our formula- 
tion for studying the electronic structure of wurzite quantum dots with cylindrical symmetry 
will be a very useful tool in this important research area. 

IV. CONCLUSIONS 

We have extended the Sercel-Vahala technique to WZ heterostructures with cylindrical 
symmetry and shown this to be an exact result. This is crucial for applications in enabling 
greatly decreased computational resources. We have studied GaN and CdSe cylindrical 
quantum dots, calculating the effect of the aspect ratio on energy levels and wavef unctions. 
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FIG. 1: The energies of the top three valence sublevels of a WZ GaN cylindrical quantum dot, 
calculated with the 3D program (circles). The solid lines are the converged values obtained with 
the 2D program. 
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FIG. 2: The energies of the lowest states as a function of the CdSe quantum dot aspect ratio. 
The linear term is not included in the Hamiltonian. Shown are the first four states with F z = 1/2 
(solid lines) and with F z = 3/2 (dotted lines). 
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FIG. 3: The energies of the lowest states as a function of the CdSe quantum dot aspect ratio. The 
linear term is included in the Hamiltonian. Shown are the first four states with F z = 1/2 (solid 
lines) and with F z = 3/2 (dotted lines). 
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